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ABSTRACT 

Simple arguments suggest that magnetic fields should be aligned tangentially to the 
surface of an AGN-blown bubble. If this is the case, charged particles from the fully 
ionised intra-cluster medium (ICM) will be prevented, ordinarily, from crossing the 
boundary by the Lorentz force. However, recent observations indicate that thermal 
material may occupy up to 50% of the volume of some bubbles. Given the effect of the 
Lorentz force, the thermal content must then be attributed to one, or a combination, 
of the following processes: i) the entrainment of thermal gas into the AGN outflow that 
inflated the bubble; ii) rapid diffusion across the magnetic field lines at the ICM/bubble 
interface; iii) magnetic reconnection events which transfer thermal material across 
the ICM/bubble boundary. Unless the AGN outflow behaves as a magnetic tower 
jet, entrainment may be significant and could explain the observed thermal content 
of bubbles. Alternatively, the cross-field diffusion coefficient required for the ICM 
to fill a typical bubble is ~ 10 16 cm 2 s _1 , which is anomalously high compared to 
predictions from turbulent diffusion models. Finally, the mass transfer rate due to 
magnetic reconnection is uncertain, but significant for plausible reconnection rates. 
We conclude that entrainment into the outflow and mass transfer due to magnetic 
reconnection events are probably the most significant sources of thermal content in 
AGN-blown bubbles. 
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1 INTRODUCTION 

Observational studies indicate that powerful radio emis- 
sion from Active Galactic Nuclei (AGN) at the centres of 
galaxy clusters is strongly related to the thermal state of 
the hot, X-ray emitting plasma that permeates the clus- 
ter. In particular, systems with short radiative cooling times 
at the cluster centre seem more likely to exhibit on-going 
star formatio n, op tical line-emission and radio emission from 



AGN (e.g. 



lBurndll99d; [Crawford et al.ll 19991; iRaffertv et~aH 
|200S| ; ICavagnolo et alJN2008l : iMittal et alj l2008h . There is 
also a significant body of theoretical evidence indicating 
that outflows from AGN play an important role in the 
evolution of their surroundings (e . g. iBinnev fc Tabod 
Benson et"aH 120031: ICroton et al.1 l200d. iBower et al " 



Short fc Thomas! 120091 ; IPuchwein et alj I2OO8I ; I Bower et atl 
20081 ; |Popei l2009). As a result, it is often assumed that AGN 
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operate as thermostats which balance and regulate the ra- 
diativ e losses of the intra-cluster medium (ICM) (e.g. iKaiseJ 
12007ft . In the standard paradigm, material cools out of the 
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ICM and is accreted by a supermassive black hole located 
in the central galaxy. This releases vast amounts of energy, 
often in the form of outflows, which couple to and heat the 
amb ient gas, thereby redu cing the cooling rate of the ICM 
(e.g. IChurazov et"ai1l2005l ). 

To understand how the energy of the outflow is actu- 
ally dissipated and heats the ICM, it is necessary to build 
an accurate model describing the interaction between the 
outflow and the ICM. In this article, we focus specifically on 
the interaction between the ICM and the bubble-like struc- 
tures inflated by AGN outflows. Observations indicate that 
these bubbles are largely d evoid of thermal material (e.g. 
iMcNamara fc Nulsenl 120071) and primarily consist of mag- 
netic fields and high-ene rgy particles (c.f. iDunn fc Fabian! 
|2004 iBirzan et all 120081 ). In order for the bubbles to ap- 
pear as depressions in the X-ray surface brightness (e.g. 
IMcNamara fc Nulsen|[2007l ). we infer that the mass of ther- 
mal material entering from the ICM must be small compared 
to the mass initially displaced during the inflation. There- 
fore, the observations seem to indicate that the bubbles do 
not mix significan t ly wi th their surroundings. In addition, 
ISanders fc Fabianl (|2007l ) argue, based on an analysis of the 
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bubbles in the Perseus cluster, that the volume fraction of 
thermal material in bubbles is no more than 50%. 

It is generally supposed that mixing is mostly inhib- 
ited by the presence of magnetic fields at the interface 
between the bu bble and the ICM. In particular, follow - 
ing the work of IChandrasekharj (|l96lh . |Pe Yound (|2003h ; 
iKaiser et aD (|2005h : [5 ursi showed that a magnetic 

field oriented tangentially to the bubble surface may pre- 
vent the growth of Rayleigh- Taylo r instabilities. This ha s 
also been sh o wn n umerically by iRobinson et alj (12004); 
lO'Neill et all (I2009T). However, s ee als o iPizzolato fc Sokeij 
( 20061 T TScannapieco fe Briiggenl (|2008h for other explana- 
tions of bubble stability. 

One possibility for the origin of such a favourable 
magnetic field orientation is that, during the expansion of 
the bubble, the ambient medium, along with its magnetic 
field, is excluded from the bubble. This results in a sheath 
around the bubble in which the field is primarily tangen- 
tial to the bubble surface, and which has a higher value 

than the ambient value, due to the effects of compression 
i it — l 

(Dc Young 2003). Nu merica l simulations of moving bu bbles 
(see [Ruszkowski et al. 2 0071 ; iDursi fc Pfrommeil [20081 have 
also shown that that magnetic field from the ICM can be- 
come draped across the front of the bubble. This ensures 
that, whatever the initial configuration of the magnetic field 
on the bubble surface, there should always be layer of mag- 
netic field that is oriented tangentially to the bubble surface. 

However, suppressing the growth of fluid instabilities 
on the bubble surface is only one aspect of the interaction 
between the ICM and the bubble. That is, the magnetic 
field on the bubble surface might provide a stable frame- 
work which prevents the bubble from fragmenting into many 
smaller bubbles, but this does not necessarily prevent ther- 
mal material entering the bubble. An important and related 
effect arises from the fact that the ICM is a fully ionised 
plasma and therefore has a high electric conductivity, while 
the AGN-blown bubble contains a magnetic field. This sce- 
nario is extremely similar to the interaction between the 
solar wind and the Earth's magnetosphere (see for e x ample 
Chapman fc Ferrard [l930l ; iKivelson fc Russelll Il995l ; iParkj 
2003). Therefore, in the idealised case, where the ICM car- 



ries no magnetic field, we expect a charge entering the bub- 
ble to experience a Lorentz force which is perpendicular to 
both its direction of motion and the magnetic field. Pro- 
vided the bubble magnetic field is oriented tangentially to 
the bubble surface, the charge will re-directed out of the 
bubble and emerge back in the ICM. The combination of 
this exclusion and the stabilising effect of magnetic fields 
creates a relatively impermeable barrier between the bubble 
and the ICM. Furthermore, the ICM is not a perfect con- 
ductor and some of the electromagnetic energy, generated 
by the interaction, will be dissipated as heat in the region 
around the bubble. In principle, this could heat the ICM 
around the bubble, or help to power the optical emission 
from cool material observed behind some bubbles. 

The plasma exclusion effect described above has signif- 
icant implications for interpreting observations of bubbles. 
In particular, any thermal material could only have entered 
the bubble by the following processes: i) thermal material en- 
trained into the outflow that originally inflated the bubble; 
ii) particle diffusion across the magnetic field on the bubble 
surface; iii) non-ideal processes such as magnetic reconnec- 



tion, which occur as a result of non-zero resistivity in the 
ICM and the interaction between magnetic fields in the ICM 
and the bubble. With reference to the firs t process, mass- 
loadi ng will decelerate the AGN outflow l|Hartauist et alj 
1986) and so may potentially explain the abs ence of strong 
shocks in the centres of g alaxy clusters (e.g. iForman et al.l 
120071 ; iGraham et"al1l2008l ). The velocity of the outflow also 
plays a role in governing the morphology of the subsequent 
structure. For example, slower outflows are more likely to in- 
flate bubble- like structures, rather than the elongated struc - 
tures produced by faster jets (e.g. ISternberg fc So ker 2008). 
Therefore, determining the origin of a bubble's thermal con- 
tent is of importance in the wider context of AGN feedback. 

The article is structured as follows: in section 2, we in- 
troduce the model which describes the interaction between 
the ICM and the bubble magnetic field. The main quantities 
derived from this model, such as the Ohmic heating rate and 
induced current are given in section 3. Section 4 provides a 
discussion of the mass transport rates associated with out- 
flow entrainment, diffusion and magnetic reconnection. We 
summarise the findings in section 5. 



2 MODEL 

In this derivation, the ICM is taken to be an unmagnetised, 
fully ionised plasma consisting only of electrons and protons. 
The bubble's magnetic field is taken to be aligned tangen- 
tially to the bubble surface, see figure 1. We also assume 
that both the ICM and bubble fluids can be described us- 
ing the equations of ideal magnetohydrodynamics. This as- 
sumption is valid in the limit that the diffusion of magnetic 
fields is much less important than their advection. The rel- 
ative importance of diffusion and advection for magnetic 
fields is often quantified using the magnetic Reynolds num- 
ber, Re m = UL/ri, where U and L are the velocity and 
length scales of interest and rj is the ele ctric resistivity. For 
a full y ionised plasma 77 is given by (e.g. iBodenheimer et al.l 
120071 ) 



rj « 260 



lnA\ / T 
~W) Vl0 8 K 



-3/2 
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where In A is the Coulomb logarithm. Therefore, 
fle m «6xl0'" ' W W 
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Since 7?e m S> 1, the diffusion of magnetic fields is unimpor- 
tant on the length scale of interest, L. More specifically, we 
can say that a magnetic field initially residing in an AGN- 
blown bubble cannot diffuse far and so is effectively frozen 
to the bubble fluid. This is the so-called 'frozen-in' approxi- 
mation, which corresponds to the flow regime known as ideal 
magnetohydrodynamics . 

The interaction between the ICM and the magnetised 
AGN-blown bubble can be understood by using a combi- 
nation of fluid and particle approaches. The fluid model 
provides a description of the conditions at the ICM/bubble 
interface, showing that in order to prevent thermal ICM 
plasma entering the bubble there must be an electric cur- 
rent flowing in the boundary. The particle description shows 
that the ICM is excluded from the majority of the bubble 
due to the Lorentz force acting on the charges, but that the 
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partial penetration and deflection of these particles is suffi- 
cient to generate the current. These methods are described 
in detail below. 



2.1 The fluid model 

Observationally, AGN-blown bubbles appear as depressions 
in th e X-ray surface brightness (e.g. iMcNamara fe Nulsenl 
meaning that they contain comparatively little mate- 
rial emitting at these energies. Recent evidenc e of bubbles in 
the Perseus cluster (ISanders fc Fa bian 2007) suggests that 
the volume fraction of thermal material is less than 50% 
meaning that the flui d pressure ac ross the boundary must 
be discontinuous (e.g. iParksl |2 003l ) . The conditions at the 
boundary must satisfy the fluid momentum equation. For a 
species, i, of charged fluid, the momentum equation is 



Pi-Qj- +Pi(Vi.V)vi 



-VP + Ji x B + piF, 



(3) 



where pi is the density of the fluid with pressure, Pi, Ji is 
the current density vector, B is the magnetic field vector 
and F is an external force per unit mass, such as gravity. 
The Ji x B term describes the reaction of the fluid to the 
magnetic field. 

Therefore, assuming that the external force, F, is locally 
unimportant, the equilibrium condition to be satisfied at the 
bubble surface is VP = Ji xB. The magni tude and dir ection 
of the induced current must then be (e.g. Parks 2003) 

Furthermore, the current calculated from equation (Q must 
also satisfy Ampere's law 

<9E 



VxB = no J + h £q 



dt 



(5) 



where fio is the vacuum permeability, and eodE/dt is known 
as the displacement current, which can be neglected. This 
must produce the J x B force needed to balance the rate of 
change of 1CM momentum. In the simplest case, where there 
is no coupling of energy or momentum across the boundary, 
these currents close on themselves. 

An estimate of the induced current density, J, requires 
a value for either VP or V x B which can be calculated 
using the particle model, described below. 



For the idealised case, the thickness of the interaction 
boundary would be approximately the plasma skin depth, 
5 = c/ujp, where c is the speed of light and lj p is the 
plasma frequency dFunaki et al.ll2007l ). However, because of 
their heavier mass, protons tend to penetrate more deeply 
into the magnetic field than electrons. This leads to charge 
separation with a negative charge layer above the positive 
charge layer. Thus, the outwardly pointing electric field re- 
strains the protons, while the electrons gain energy (e.g. 
iBaumiohann fc Treumannll 19961 ). 

The thickness of the boundary is therefore considered 
to be larger than <5, being roughly the proton gyration radius 

r - ^£ (6) 

where m p is the proton mass, v p is the velocity perpendicular 
to the field line, e is the electric charge, and B is the mag- 
netic field strength. Since the bubble ascends subsonically, 
the particle velocity is largely governed by its thermal mo- 
tion. Then, for a thermal proton, the mean velocity is given 
by v p = (Sfer/rap) 1 ' 2 , where ks is the Boltzmann con- 
stant and T is the temperature. Therefore, assuming that 
the electrons and protons have the same temperature, the 
ratio of the thermal proton and electron gyration radii is 
(m p /m c ) 1/2 . 

The total current density in the y-direction is the sum 
of the individual currents generated by the protons and elec- 
trons: |J| = Ji = (n p ev p + n c ev a ), where n p and n c are 

i— p,e 

the number density of protons and electrons, respectively. 
v p and v e are the proton and electron velocities. 

Then, using Ampere's law, the magnitude of the in- 
duced current density is roughly 

£ W « £ (7) 

where r g ,i is the gyration radius of the ith particle species 
and P is the pressure of the ith species. 

We also note that the current flows so as to exclude the 
bubble's magnetic field from the ICM, which increases the 
field strength just inside the boundary region. The numeri- 
cal value of the field inside the boundary is then twice the 
initial value. As a result, the current layer confines the bub- 
ble's magnetic field to the region inside the boundary and 
so partly governs the size of the bubble. 



2.2 Particle model 

In the scenario depicted in figure 1, the incident ICM par- 
ticles move in the negative a;-direction with a velocity of 
magnitude u. The magnetic field points in the ^-direction, 
with a magnitude B so that the Lorentz force acting on 
an incident proton is in the y-direction, with a magnitude 
|F| = euB. Once moving in the y-direction, the Lorentz force 
on the proton is then directed in the a;-direction causing 
it to be excluded from the bubble. The transverse velocity 
component (in the y-direction) that arises during the half- 
orbit generates an electric current in the boundary which 
prevents the magnetic field from entering the ICM, and is 
usually referred to as the Chapman- Ferraro current (e.g. 
IChapman fc Ferrarol Il930l ; iFunaki et all 120071 ). Therefore, 
once back in the ICM, a charge experiences no further de- 
flection. 



3 RESULTS 

Since the electric conductivity of the ICM is finite, some 
of the electromagnetic energy generated by the interaction 
will be dissipated as heat in the region around the bubble. In 
this section, we present a calculation of the expected heat- 
ing rate. The magnitude of the induced current is also of 
interest, primarily for comparison with th e values required 
in m odels of current-dominated jets (e.g. iNakamura et al.l 
2006). Accordingly, a calculation for the total current flow- 
ing across the surface of a bubble is also given in this section. 

3.1 Heating rate 

From equations ([6]) and @, it is clear that particles with 
smaller gyration radii contribute a larger current density. 
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Figure 1. Schematic showing the interaction of the ICM with the magnetic AGN- blown bubble, based on Baumjohann & Treumann 
(1996) and Parks (2003). 



Because of this effect, the induced current density due to 
electrons is a factor of {m p /m B Y^ 2 greater than for protons 
and has a magnitude 
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Although the electric conductivity of the ICM is extremely 
high, it is not infinite. Consequently, the induced current 
dissipates electromagnetic energy at a rate e = r)\ J| 2 , where 
rj is the electric resistivity given by equation |T]). This process 
is known as Joule or Ohmic dissipation. Due to the J term, 
the Ohmic dissipation rate is a factor (m p /m e ) higher for 
electrons than protons and can be written 



■ri\J* 
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erg s cm 



(9) 

A plasma with a temperature of T ~ 10 8 K and a number 
density n ~ 0.01 cm" 3 primarily radiates energy by ther- 
mal bremsstralung at a rate n 2 A(T) ~ 10 -27 erg s _1 cm -3 . 
As can be seen, from equation ((9J the Ohmic dissipation 
rate within the boundary region is a factor of ~ 10 9 greater 
than the radiative losses due to thermal bremsstrahlung in 
the ICM surrounding the bubble. This is surprisingly large, 
suggesting that Ohmic dissipation might have an observable 
effect. To determine whether this is true it is necessary to 
calculate the total rate at which energy is dissipated. 

Assuming that the bubble is spherical, the volume 
within which the electromagnetic energy is dissipated is 
(4/3)7r[(rb +r g ) 3 — rjj] ~ 47rr 2 r g . Because of the extra factor 
of r g , the volume-integrated dissipation rate for electrons is 
a factor (nip/me) 1 ' 2 greater than for protons, so electron 



dissipation still dominates the total Ohmic dissipation rate 

2 / r> \ 3 



E « 47rr b Y g , c e «2xl0 
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(10) 



It is also instructive to demonstrate how the heating rate 
may vary as the bubble evolves. In a stratified atmosphere, 
the rising, buoyant bubble expands to maintain pressure 
equilibrium with its surroundings. For the present case, 
pressure equilibrium requires that B 2 /(2fio) — P, where 
P is the ambient pressure. Then, assuming the expansion 
is adiabatic, the bubble volume, Vb, is related to the am- 
bient pressure, P, by PV£ = F\>,oV^ , where Pb,o is the 
ambient pressure where the bubble was inflated, Vb,o the 
initial bubble volume and F is the adiabatic index of the 
bubble. A magnetically-dominated bubble can be described 
by a relativistic equation of state, for which T — 4/3. Fi- 
nally, the ICM pressure can be approximated by a /3-model: 
P = Po[l+ (z/zo) 2 ] ' 3 , where Po is the central pressure, z is 
the distance from the cluster centre, zo is the scale height of 
the distribution, and typically /3 « 3/4. Assuming the bub- 
ble was initially inflated at the cluster centre, and substitut- 
ing for both rh and B, the heating rate can be re-written 
as 

3/2 



E « 2 x 10° 
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The thickness of a shell around the bubble that will be 
heated by Ohmic dissipation can be estimated from E c — 
A(T), where n 2 A(T) are the radiative losses due 
to thermal bremsstrahlung. Using similar quantities to those 



above, it can be shown that Ar 
will not be observable. 
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Therefore, Ohmic dissipation should be insignificant 
compared to the total radiative ICM cooling rate due to 
thermal bremsstrahlung and is unlikely to have a noticable 
effect even in a thin shell around the bubble. Furthermore, 
even if the conductivity was drastically reduced, the small 
volume occupied by the induced current means that the to- 
tal heating rate will always be negligible compared to the 
radiative losses of the ICM. 



3.2 Total current 

The current generated by the Lorentz force flows in a thin 
layer across the surface of the bubble. Assuming, again, that 
the bubble is spherical with a radius rt,, the total available 
area over which the current can flow is 47r[(Yb + r g ) 2 — r^] « 
8nrbr g . Then, using equation (O, the total current flowing 
around a typical bubble is 



/ « 87T7"b Jir s i = ^^-rtB 
^ Ho 



6X10 1 



rh 



B 



K 5kpc J 

(12) 

and is independent of the gyration radius. 

The magnitude of the current given in equation (|12[) is 
somewhat uncertain because although the current loop prob- 
ably closes somewhere at the rear of the bubble, it is impos- 
sible to say precisely where. Because of this uncertainty, the 
current may only flow over a fraction of the bubble surface, 
thereby reducing the total current. However, the reduction 
is unlikely to be more than an order of magnitude. 

Nevertheless, the value of the current in equation (|12l) is 
almost id entical to the values expected in current-d ominated 
jets (e.g. iNakamura et al.|[200a iDiehl et al.ll2008l ). In fact, 
this is not a coincidence - the electric current and bubble ra- 
dius associated with a magnetic field in pressure equilibrium 
with its surroundings are related by Ampere's law. How- 
ever, the values arise from different starting points. For th e 
current-dominated jet, described bv lNakamura et"al"1 (2006), 
the constant current generates a magnetic field; the bubble 
radius can then be calculated assuming the field is in pres- 
sure equilibrium with the ICM. In contrast, we have cal- 
culated the induced current for an assumed magnetic field 
(in pressure equilibrium with its surrundings) and bubble 
radius. The main difference is that the current estimated in 
this way is not constant, and decreases as the bubble ascends 
through the ICM. 



4 DISCUSSION 

In the idealised model of the interaction between the ICM 
and the AGN-blown bubble, described above, there is no 
mass or energy transfer across the interaction boundary (as 
well no drag force exerted on the bubble). If this is true, then 
it is not clear there should be any thermal material in the 
bubble at all, while recent observations indicate that thermal 
material may occupy up to 50% of the volume of som e bub- 
bles in the Perseus cluster l|Sanders fc Fabian|[2007h . There 
are three main possibilities which can explain this discrep- 
ancy: 

1) Thermal material from the ICM was entrained into 
the outflow that inflated the bubble. It may also be impor- 
tant to note that if a supersonic outflow entrains material, it 



A, 



will b e decelerated so that its Mach number tends towards 
unity l|Hartquist et al.lll98rj ). This effect may explain preva- 
lence of_shocksjritherCMwi^ close to unity 
(e.g. iForman et ai1l2007l : iGraham et al.| [2008) . as well as the 
origin of thermal material in the bubble. 

2) The second possibility is that charged particles dif- 
fuse into the bubble. Due to the initial expansion and the 
effect of magnetic draping, the magnetic fields will be pri- 
marily tangential to the bubble surface. Therefore, in order 
to enter the bubble, particles must diffuse across the field 
lines. Unlike diffusion parallel to the magnetic field lines, 
which can occur rapidly, diffusion across field lines is signif- 
icant ly reduced by many orders of magnitude (c.f. lEnfilinl 
2003). In fact, the di ffusion coeffici ents perpendicular to the 
field lines derived by lEnfilinl l|2003h are more than 10 orders 
of magnitu de less than those empl o yed in numerical sim- 
ulatio n s bvlRuszkowski et al.l (|2007l ); iMathews fc Brighentil 
l|2008h ; lMathewsl (120091 ). We also note that, since the bub- 
ble ascends at a velocity, w, with respect to the material 
entering it, the transfer of mass is accompanied by a retard- 
ing force, —Mw, where M is the rate at which mass enters 
the bubble. Depending on the magnitude of the mass flux, 
the retarding force may play an important role in govern- 
ing the terminal velocity of the bubble. The magnitude of 
the diffusive mass flux, therefore, has potentially important 
consequences for interpreting observations of bubbles. 

3) Magnetic reconnection occurs in a medium with non- 
zero resistivity (e.g. the ICM) when oppo sitely directed mag- 
netic fields are brought close together (|Bodenheimer et al.l 
120071 ). Importantly, when the magnetic fields of the ICM and 
bubble become interconnected, the field attains a finite com- 
ponent, B n , normal to the bubble surface, allo wing plasma 
to flo w across the bubble/ICM boundary (e.g. iPaschmannl 
Il997l) . Transfering mass in this way would deliver a retarda- 
tion to the bubble, and may also generate convective motions 
within the bubble that could affect its radio emission. 

Regarding point 1), material may be entrained into an 
outflow due to the growth of Kelvin-Helmholtz (K-H) insta- 
bilities at the boundary betw een the outflow and ambient 
medium (e.g. iDe Yo ung 2006). These instabilities form as 
a result of the velocity shear between the two fluids and 
have b een shown t o occ ur in magnetohydrodynamic (MHD) 
flows (|Rvu et all |2000D as well as purely hydrodynamic 
flows. If the masse s entrained by hydrodynamic outflows , 
10 7 - 1O 9 M (e.g. IDe Young||l986l ; [Reynolds et al.ll2002l ). 
are applicable to real AGN outflows it is possible that this 
mechanism can explain the the rmal content of AGN-b lown 
bubbles. However, recent work (|Nakamura et all 2007) sug- 
gests that magnetic tower jets do not develop K-H insta- 
bilities. Thus, if AGN outflows follow this description, the 
thermal content of bubbles must arise from either cross-field 
diffusion or magnetic reconnection, discussed below. 

The diffusive mass flux, per unit area, can be written 
F — DV p, where D is the diffusion coefficient and Vp is 
the density gradient. The total mass flux into the bubble 
is then obtained by integrating F over the surface of the 
bubble, which, for a spherical bubble gives MdiB = ^r\F . 
Furthermore, assuming that the bubble is initially devoid 
of thermal content and the bubble/ICM boundary has a 
thickness of roughly one proton gyration radius, the density 
gradient, V p, can be no greater than p/r gtP . As a result, the 
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rate at which mass diffuses into a spherical bubble cannot 
exceed 4ivr^Dp/r gtP . 

A convenient way of expressing the importance of dif- 
fusion is the 'refilling' timescale over which the diffusive in- 
flux of material becomes comparable to the mass of ICM 
material displaced by the bubble, r^je = Miis/iWdiff . Us- 
ing Mdis = (4/3)7rrjjp therefore, gives an upper limit on the 
refilling timescale due to diffusion 



rdifi ^5^^ al (l0«cm2s-i) (rffe) (13) 

X (l^c) (wk) Gyr ' 



Equation (|13|l shows that a typical cross-field diffusion coeffi- 
cient of D ~ 10 16 cm 2 s _1 is required to transfer a significant 
quantity of ICM material over a bubble lifetime of 100 Myr. 

The diffusion coefficient must fall within the range 
^cias ^ D ^ Db, where D c i as and Db are the classical an d 
Bohm diffusion coefficients, respectively (|lkhsanovll200lT) . 
Classically, the diffusion of charges across magne tic field 
lines arises from non-zero plasma resistivity (e.g. ISpitzerl 
1 19621 ; iGoldston fc Rutherford! 1 1995f ). that is, it scales with 
the collision rate of electro ns with protons, v e ,p, and i s 
roughly D c i as w u e , p r^ <e (e.g. IGoldston fc Rutherfordlll995h . 
where r giC is the electron gyration radius. For typical clus- 
ter values, -D c i as ~ 10cm 2 s _1 and, therefore, will not lead 
to significant mass transfer from the ICM to the bubble on 
a realistic timescale. 

The maximum diffusive mass transfer rate occurs if dif- 
fusion is gove rned by drift-dissipative ins tabilities - Bohm 
diffusion (e.g. llkhsano^l200ll ; lEnfilirJl2003l) . for which D B « 
k B T/(16eB) w 5 x 10 15 (T/10 8 K)(S/10/x G) _1 cm 2 s" 1 . 
Therefore, mass transfer due to Bohm diffusion would be 
able fill a bubble on a timescale of 100 Myr. Despite be- 
ing an extreme upper limit, the Bohm limit does, neverthe- 
less, provide an estimate of the minimum possible timescale 
a bubble can survive before diffusion-driven mass transfer 
completely refills it. 

Perhaps more realistically. lEnfilinl (|2003h calculated dif- 
fusion coefficients for cosmic rays crossing a turbulent mag- 
netic field to be ~ 10 15 cm 2 s" 1 , using B = 10 fiG. For ther- 
mal particles with an energy of ~ lOkeV, the diffusion co- 
efficient is reduced to D ~ 10 11 cm 2 s~ 1 . In this case, the 
time taken to completely refill the bubble would be 10 12 yr. 
Therefore, unless diffusion proceeds at the Bohm rate, it is 
unlikely to contribute a significant propertion of the thermal 
content of bubbles. 

If magnetic reconnection occurs, the plasma transfer 
rate is determined by the dimensionless reconnection rate 
commonly defined as £ = B n /B. However, even for the 
Earth's magnetosphere it is difficult to accu rately determine 
B n , though typical values indicate C ~ 0-1 l|Elsner fc Lambl 
1 19841 ; |Paschmannlll997l ; IlkhsanovllioOll l. The mass flux en- 
tering the bubble is then given by M Tec = irr^pwC, wh ere w 
is the ascent velocity of the bubble l|Paschmann| [l997). The 
timescale for refilling a bubble as a consequence of magnetic 



reconnection is then 

_ M dis 4 r b / r b \f w y 1 

TlCC ~ Af rcc ~ 3K VSkpcAlO^cms-iy' (14j 

According to equation (|14[) , it seems physically possible 
that magnetic reconnecti on could explain the observations of 
ISanders fc Fabianl (|2007h . Furthermore, the transfer of mass 
would also be accompanied by heating as the energy density 
of the magnetic field is transferred to the plasma. However, 
the significance of this channel of mass transfer depends cru- 
cially on the reconnection rate, £, which is highly uncertain. 



5 SUMMARY 

The aim of this article has been to describe one of the main 
processes which prevents thermal material from the ICM 
entering AGN-blown bubbles. As a direct consequence, it is 
possible to highlight the mechanisms which must occur if a 
bubble is to display thermal content. The main findings are 
summarised below: 

i) The Lorentz force causes a charged particle entering 
the bubble to execute a half-orbit so that it is re-directed 
back into the ICM and excluded from the bubble. The trans- 
verse motions of the deflected charges generate electric cur- 
rents along the bubble surface which prevent the magnetic 
field from entering the ICM (e.g. IChapman fc Ferrardll930l : 
iFunaki et all 120071 ). Therefore, once back in the ICM, a 
charge experiences no further deflection. The charges pene- 
trate roughly one proton gyration radius before being di- 
rected back into the ICM. This exclusion is different to 
the stabilising effect of magnetic fields (e.g. IChandrasekharl 
Il96ll ; |Pe Youndl2003l ; iKaiser et af1l2005l ; |Pursill2007l ) which, 
strictly speaking, only explains why the bubbles do not frag- 
ment into smaller structures due to Rayleigh- Taylor insta- 
bilities. The combination of both effects provides an effective 
surface tension which, in the idealised model, prevents ma- 
terial from the ICM entering the bubble. 

ii) Ohmic dissipation of the induced current heats the 
ICM adjacent to the bubble at a rate ~ 10 35 - 10 37 ergs _1 ; 
this is insignificant compared to the radiative losses of the 
ICM. The total induced current on the bubble surface has a 
typical magnitude of I ~ 10 18 A, which is similar to that 
required to produce bubbles from curr ent-dominated jets 
jDiehl et al.ll2008l; INakamura et ai1l2006l ). 

iii) ISanders fc Fabianl (|2007l ) recently presented obser- 
vations indicating that thermal material may occupy up to 
50% of the volume of some bubbles in the Perseus cluster. 
Given the effect of the Lorentz force, the most likely expla- 
nations for the observed thermal content are that: mass is 
entrained into the outflow that inflated the bubble, charged 
particles diffuse rapidly across the bubble/ICM boundary, 
or magnetic reconnection on the bubble surface allows ma- 
terial to enter the bubble. We find that, unless the AGN 
outflow behaves as a magnetic tower jet, entrainment into 
the outflow ca n supply sufficient t herma l material to explain 
the results of ISanders fc Fabianl |2007). However, for mag- 
netic tower jets, th e outflow entrainment rate is expected 
to be insignificant (INakamura et al1l2007h . The cross- field 
diffusion coefficient required for thermal material to fill a 
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typical bubble is ~ 10 16 cm 2 s _1 . This value is anomalously 
hig h compared t o predictions from turbulent diffusion mod- 
els (|En6rinll2003h . but is comparable to the Bohm diffusion 
coefficient. Given that the Bohm limit is an extreme upper 
limit, it seems unlikely that cross-field diffusion contributes 
significantly to a bubble's thermal content. Finally, the mass 
transfer rate due to magnetic reconnection is somewhat un- 
certain, but significant for plausible reconnection rates. We 
conclude, based on current models, that entrainment into 
the outflow and mass transfer due to magnetic reconnection 
are probably the most significant sources of thermal content 
in AGN-blown bubbles. 
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